Critical properties of the S—l spin dimer compound Ba 3 Mn 2 08 
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BasMn2 08 is a hexagonally coordinated Mn 5+ 5=1 spin dimer system with small uniaxial 
single-ion anisotropy. 135,137 Ba NMR spectroscopy is used to establish the lower critical field 
H c i of distinct field-induced phases for H || c,H _L c, and measure the longitudinal (Mi) and 
transverse (Mt) magnetizations in the vicinity of the quantum critical point (QCP). Mm (T, H c i), 
Mi±(T, Hd) are reproduced by solving a low-energy model for a dilute gas of interacting bosons. 
M t ||(T -> 0,H = H cl ) (M t± (T -> 0,H = H cl )) follows the expectation for a BEC (Ising-like) QCP. 

PACS nos. 75.45.+j,75.40.Cx,75.40.-s,76.60.Cq 



Recent investigations of field-induced phases in £=1/2 
magnetic insulators typify the opportunities for studying 
the problem of Bose Einstein condensates (BEC's) specif- 
ically [H, and quantum criticality more generally. In 
spin-dimcr, and other spin-gapped systems, the ground 
state is a singlet while the lowest excited states are 
a mode of triplet excitations @, 0]- The magnetic 
field tunes the chemical potential for triplet excitations 
through zero at the critical field H c \ producing a con- 
trolled density of triplets, that can either condense or 
crystallize into a superlattice depending on the balance 
between kinetic and potential energies 0, [B[ . The Hamil- 
tonian has U(l) rotational symmetry in the idealized 
case, and this symmetry is spontaneously broken in the 
condensed phase with the development of a finite trans- 
verse magnetization M t . 

From what is known about the spin-dimer system 
Ba3Mn208 [H, these conditions hold for H || c Q. How- 
ever, the evolution of the phases in a magnetic field is 
known to deviate from the simplest 5=1/2 isotropic case 
in a number of ways Two magnetization plateaus 

with (S Z }=1 (per dimer) and (S z )=2 are observed as a 
result of the 5=1 state of the Mn 5+ ions @, G3- In 
addition, a small single-ion uniaxial anisotropy is under- 
stood to produce new boundaries in the ordered phases 
for H tilted from the -axis. While this anisotropy is not 
relevant for H || c, its influence is most prominent for 
H _L c, where there is evidence for an additional phase 
II, stabilized only near H c ±± and the other three critical 
fields. Further, the hexagonal coordination of the layers 
leads to geometric frustration. The near-neighbor trans- 
verse spin components would be rotated by cv=120° in 
an isolated triangular layer ll|, ll2j. However, interlaycr 
coupling (Fig. Q) leads to a -> 120° + e with e - 9°), 
because incommensurate spin ordering partially releases 
the interlayer frustration. 

Presented here are results of 135 ' 137 Ba NMR spec- 




FIG. 1: (color online) a) Ionic arrangement of Ba3Mn20s 
showing location of Ba sites relative to the Mn ions; oxygen 
ions are not shown. The exchange couplings are indicated by 
solid (Jo) and dashed (Ji, J 2) lines, b) Illustration of the two 
chiralities, with the transverse spin directions denoted by the 
arrows, with 01 = 129° (see text). 



troscopy studies in the high symmetry phase near H c i, 
and in the ordered phases I (H || c, H _L c) and II 
(H _L c) . The NMR shifts are used to establish the mag- 
netization as a function of temperature at H =H C \. For 
the case H || c, the result is consistent with the expecta- 
tions for a BEC-QCP, i.e., M(T -> 0, H cl ) ~ T 1 ' 5 . Both 
the universal and the non-universal (T >100mK) regimes 
are quantitatively described using an effective low-energy 
theory for a dilute gas of bosons. Quantitative differ- 
ences are observed for H _L c, in agreement with the 
expectation for the Ising-like (Z2) broken symmetry of 
phase II 0. Further, we investigate the longitudinal and 
transverse magnetization (Mi, Mt) of the ordered phases, 
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and establish that the field dependence of the transverse 
magnetization follows the expected mean-field behavior 
(D = d + z > 4). We also conclude that the line of 
transitions dividing I /II is discontinuous. 

The maximum temperature of the ordered phases 
is T m =0.9K Q, so the measurements reported here 
were performed on a single crystal placed inside the 
mixing chamber of a dilution refrigerator. 135 > 137 Ba 
( 135 ' 137 I=3/2) NMR spectroscopy was performed in mag- 
netic fields H <120kOe using a top-tuned configuration. 
The platform holding the sample and coil is rotated by an 
Attocube piezoelectric motor. At the higher fields avail- 
able at the NHMFL, we used a bottom-tuned 3 He system. 
The diagonal hypcrfinc couplings were determined by 
comparing high temperature measurements of the shift 
(T >20K) to susceptibility measurements @. Orbital 
and quadrupolar couplings were determined from the 
shifts measured at the lowest temperatures for H < H c \ . 
Some of the NMR parameters are summarized in Table 
III 

The £2 = 1 triplet excitations become gapless at H c \, 
while the other two triplets and the quintets have a gap 
of order Jo. Therefore, only S z — 1 triplets can be 
created at low energies for H ~ H c i. These triplets 
are hard-core bosons with an effective chemical poten- 
tial p = g vv pB{H — H c i) (y = {a, 6, c}). The effective 
Hamiltonian for H || c (y = c) or H _L c [y = a, b) is: 



q q 

2N 



q,k,k' 



(1) 



where e q = ^Af + fAiJq - ^/A 2 + f Ai J Q , Q is the 
wave- vector that minimizes e q , a q = gJ7(Q), Ai = 
1.65meV is the single-dimcr singlet to S* = 1 triplet gap, 
J q = 2( h - Jshl + f 7 q + f 7 q , and 



7 q = cos q 3 + cos (q 3 -q x )+ cos (q 3 -q 2 ), 
= cosq! + cos g 2 + cos (qi — q 2 ), 



7 q = cos (q 3 - q 2 + qi) + cos (q 3 - qi + q 2 ) 



cos (q 3 - qi- q 2 ). 



(2) 



The relative magnitude of the exchange anisotropy g = 
0.088 was obtained by fitting the difference between the 
values of H cl u and H cl± : AH cl = 89.3 - 86.4 = 2.9kOc. 





A-aa 




137 i^q (MHz) 


Ba(l) 


0.26(1) 


0.35(2) 


54.7(2) 


Ba(2) 


0.18(1) 


0.11(1) 


10.8(2) 



TABLE I: Selected NMR parameters for the two sites shown 
in Fig. □ v Q = e 2 qQ, with 135 Q ( 137 Q)=0.18 (0.28) 
xlO -24 cm 2 . The hyperfine coupling constants are reported 
in ^ S /Mn 5+ . 



The exchange constants [l3( and the g-factors are Ji = 
0.118mcV, J 2 - J 3 = 0.114meV, J 4 = 0.037meV, g cc = 
1.98 and g aa = 1.97. The effective repulsive interaction 
v o = ToCQjQ) results for summing the ladder diagrams 



for the bare interaction vertex Vq [14 1 



8tt 3 Kq ~ p 



r P (k,kQ 

e k+p 

(J2 + J3K 



(3) 



r q (k,k') = v q - 

For Ba 3 Mn 2 8 , we have V^ = U + ( J 2 + J 3 )7 q + 
irTq' wnerc V ~ * 00 comes from the hard-core repulsion, 
while the rest of the terms correspond to the off-site re- 
pulsive interactions the result from the Ising terms of the 
inter-dimer exchange couplings. By solving Eq.Q, we 
obtain v = 0.9meV for J 2 + J 3 = 2.82K. The value of 
J2 + J3 = is obtained by fitting H c2 ~ 27T for H || c. 
We note that the second term of H breaks the U(l) sym- 
metry associated to the conservation of the total number 
of bosons (A/;) for H _L c. This term is originated by 
the effective exchange anisotropy found in Ref. Q ■ Con- 
sequently, we expect an Ising-like (broken Z 2 ) quantum 
phase transition (QPT) for H 1 c in contrast to the 
BEC-QPT that occurs for H || c. 

Magnetization results for H || c are shown in Fig. 
[2^ for several fields near H c \. The curve measured at 
H = H c i = 89.3kOe is consistent with the expectation 
M^||~ T 3 / 2 for T -> 0. The red line is the result of a 
Hartree-Fock decoupling of the last term of Ti. that in 
the disordered phase has the effect of renormalizing the 
chemical potential /x e ff = ji — 2vop (p is the density of 
bosons) [2j. Calculations for field values differing from 
Hci, H = H cX - 1.3kOe (blue) and H = H cX + 1.5kOe 
(green), also match the NMR shift data well. There is 
a 20% disagreement if only the hard-core repulsion is in- 
cluded in Eq. ©. 

When the applied field is rotated to the ab plane, the 
ordered phase II bordering the paramagnet is believed to 
be Ising-like, with transverse spins confined to the c di- 
rection. In this case, we expect M^j_~ T 2 for T — > and 
H = H c ix- The NMR results in Fig. [^b agree well with 
the mean field treatment of Ti. The anisotropy term also 
has the effect of lowering the critical field H cX . In con- 
fining the transverse spins to the c-axis, the energy gain 
associated with the broken symmetry is reduced slightly, 
and consequently Mi±<M^. The outcome is consistent 
with the anisotropy parameter D=32jxeV as established 
by electron paramagnetic resonance 15 1. 

Turning to the transverse magnetization, Fig. [3] shows 
two field-swept spectra in the condensed phase for H || c. 
The blue segment in the inset shows the location in the 
phase diagram where these spectra were recorded. The 
spectrum of the Ba(l) site is well reproduced by assum- 
ing a simple plane wave incommensurate modulation of 
the longitudinal field (dashed red line). The functional 
form is independent of H for H || c; only the spectral 
width and overall shift vary. The Ba(2) spectra are com- 
posed of two parts of equal intensity, and there appears 
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FIG. 2: (color online) (a) M;|j vs. T for selected magnetic 
fields close to i3" cl ii . The inset shows a sequence of spectra 
recorded at different temperatures. The solid and dashed 
curves are from ladder diagram calculations (see text), (b) 
The same as (a), for H _L c. 

to be an asymmetry in the line shape. The asymmetry 
is reduced for fields close to H c \ (inset of Fig. so 
the lineshapes are modelled by plane wave incommensu- 
rate structures with different widths and shifts. A small 
systematic error is introduced by varying magnetic field 
rather than frequency, so all subsequently shown spectra 
were obtained by sweeping the frequency at constant H . 

The incommensurate modulation arises from the in- 
tcrlaycr coupling and the off-diagonal hyperfine coupling 
makes its detection possible. In comparing the broad- 
ening of central and satellite transitions, the magnetic 
origin of the broadening is confirmed [l6(. The two in- 
equivalent sites for Ba(2) are associated with triangu- 
lar coordination of the Mn sites with opposing chiral- 
ity (±) of the transverse order. In the Ba(l) case, the 
point group symmetry is (C^): each Ba(l) site is situ- 
ated equidistant from three Mn sites in each of two dimer 
layers. The two Mn triangles have opposite chirality in 
the condensed phase, and by symmetry there is no dis- 
tinction of the Ba(l) sites except for the phase modula- 
tion resulting from the incommensurability. The Ba(2) 
sites have lower (C3) symmetry, with the three nearest- 
neighbor Mn forming a triangle of specific chirality. 

The spectrum broadens and shifts near H c x- FigHk 
shows the frequency-swept Ba(2) spectra for the central 
(1/2^-1/2) transition, collected for a sequence of fields 
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FIG. 3: (color online) Field-swept spectra of 137 Ba(2) satel- 
lite and 137 Ba(l) central transition. The dashed line is the 
expectation for a plane wave incommensurate hyperfine field 
modulation. The partial T/H phase diagram in the inset re- 
sults from NMR relaxation data (solid red diamonds) UM 
and thermodynamic measurements (open squares) in Ref. 
The blue line marks the temperature and field range for the 
main panel. 



close to H c \. The field dependence of M t (H) is generated 
by using the model of an incommensurate modulation of 
the hyperfine field (see Fig. [4p). The results obtained 
from both Ba(I) and Ba(2) sites are included in the fig- 
ure. The solid and dashed lines are guides based on the 
approximate empirical relationships, Mi « 9vb(h-h c1 ) 

c2 "cl 

M t = [Miil-Mi)} 1 / 2 . The first expression linearly inter- 
polates between the critical fields, while the second is the 
simplest (two level approach) mean field result expected 
for the condensed phase I. We note that the mean field 
exponents, M/ cx (H - H cl ) and M t cx (H - iJ c i) 1/2 , are 
correct because the effective dimension D = d + z = 5 
exceeds the upper critical dimension. 

Fig [5] shows the Ba(2) spectra near H c i±. Unlike for 
H || c, the spectra are distinctly asymmetric at all fields 
in the range of phase II. Still, the spectrum remains rel- 
atively simple up to H wI00-I05kOc. The spectrum be- 
comes particularly complicated over a range of fields ex- 
tending to approximately II5kOe, then simplifies once 
again, for H > 115kOe. We take this observation as ev- 
idence for a line of first order phase transitions dividing 
phases II and I, consistent with recent neutron scattering 
experiments carried out independent of this work [f7j . 

It is useful to consider the prediction for the nature of 
phase II to interpret the significance of the lineshapes. 
At low temperatures and H > H&, the transverse mag- 
netization oscillates with wavevector Q out of the plane, 
M t || c. The longitudinal component is modulated with 
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FIG. 4: (color online) a) Frequency-swept spectra for Ba(2) 
site near to H c \. The inset shows a comparison of the result 
to two independent plane wave incommensurate modulations 
of the hyperfine field. The two modulations correspond to 
Ba(2) locations of positive (+) and negative (-) chiralities. b) 
Longitudinal and transverse magnetization inferred from the 
spectra. 
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FIG. 5: (color online) Ba(l) spectra for H _L c close to 
-Hcix=86.5kOe. Lineshape changes in the range of fields 
102kOe< H <115kOe are consistent with a line of first or- 
der phase transitions dividing phase I and II. 



wavevector 2Q. The behavior is parameterized as 



S^ = {-lY-j=Bm26coB{v r Ql S%=Q (4) 

1 1 
- ± - 

4 4 



= - A ± ^\Jcos 2 29 + sin 2 20sin 2 (r J • Q), (5) 



where j is a dimer index, \i = 1, 2 labels the two sites of 
each dimer, H _L cdefincs the z-axis, and 9 = [0, it]. The 
results of model simulations indicate that the lincshapes 
for H <105kOe are qualitatively consistent with phase 
II, provided that both isotropic and anisotropic hyperfine 
coupling terms are included [l6| . 

The 135 ' 137 Ba spectroscopy reported here summarized 
the behavior around the critical point H = H c \ for two 
directions of H. For the longitudinal magnetization, 
the data is well described by including interdimer (near- 
neighbor) repulsions in the ladder calculation. In particu- 
lar, the BEC universality class applies for the case H J| c, 
whereas Ising-like criticality applies to H _L c. Our data 
also establishes the presence of incommensurate order pa- 
rameters. However, some parameters, which are relevant 
to the analysis of the Ba(2) site hyperfine coupling for 
H || c, remain unknown. For example, the shifts of the 
sites located near Mn triangles of opposite chirality are 
different in the condensed phase, and it is not clear why 
this should be the case. A possible explanation could 
originate with spin-orbit interactions indirectly impact- 
ing the hyperfine fields in Ba(2) sites of different chirality. 
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